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ABSTRACT 

We present new detections of the CO(5-4), CO(7-6), [Cl]{^ Pi-^ P(i) and [Gl]{^ P2-^ Pi) 
molecular and atomic line transitions towards the unlensed, obscured quasar AMS12 
(z = 2.7672), observed with the Institut de Radioastronomie Millimetrique (IRAM) 
Plateau de Bure Interferometer (PdBI). This is the first unlensed, high redshift source 
to have both atomic carbon ([CI]) transitions detected. Continuum measurements 
between 70 /xm and 3 mm are used to constrain the far infrared (FIR) spectral energy 
distribution (SED), and we find a best fit FIR luminosity of logio[LFiR/L0] = 13.5 ± 
0.1, dust temperature Td = 88 ± 8 K and emissivity index /3 = 0.6 ± 0.1. The highly- 
excited molecular gas probed by CO(3-2), (5-4) and (7-6), is modelled with large 
velocity gradient (LVG) models. The gas kinetic temperature Tq, density n(H2), and 
the characteristic size rg, are determined using the dust temperature from the FIR 
SED as a prior for the gas temperature. The best fitting parameters are Tq = 90 ±8 
K, n(H2) = 103-9±o.i (,j^-3 g^j^j ^ 0.8 ± 0.04 kpc. The ratio of the [CI] lines gives a 
[CI] excitation temperature of 43 ± 10 K, indicating the [CI] and the high-excitation 
CO are not in thermal equilibrium. The [CI] excitation temperature is below that 
of the dust temperature and the gas kinetic temperature of the high-excitation CO, 
perhaps because [CI] lies at a larger radius where there may also be a large reservoir 
of CO at a cooler temperature, perhaps detectable through the CO(l-O). Using the 
[CI] i^Pi Po) line we can estimate the strength of the CO(l-O) line and hence the gas 
mass. This suggests that a significant fraction (^ 30%) of the molecular gas is missed 
from the high-excitation line analysis, giving a gas mass higher than that inferred from 
the assumption that the high-excitation gas is a good tracer of the low-excitation gas. 
The stellar mass was estimated from the mid-/near-infrared SED to be Mj, ~ 3 x 10^^ 
Mq . The Eddington limited black hole mass is found from the bolometric luminosity to 
be M. ^ 1.5 X 10^ Mq. These give a black hole - bulge mass ratio of M,/M* ^ 0.005. 
This is in agreement with studies on the evolution of the M,/M^ relationship at 
high redshifts, which find a departure from the local value ~ 0.002. We discuss the 
implications for the evolution of the black hole in AMS12 and its host galaxy. 

Key words: galaxies: active - galaxies: high redshift - quasars: emission lines - 
quasars: individual AMS12 
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1 INTRODUCTION 

Studies into the gas and dust in high redshift (z ^ 2) quasar 
hosts allow us to observe an important epoch in galaxy for- 
mation. Determining the physical properties of the gas and 
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dust allow for the characterisation of these galaxies, provid- 
ing comparisons to the nearby Universe. 

The gas and dust are thought to be heated by nearby 
star formation (SF) (i.e. heating by young OB stars), and 
possibly by the AGN itself. Investigations into possible AGN 
heating of the gas and dust are important to distinguish the 
SF and AGN contributions to the far-infrared luminosity 
(ipin), and the CO luminosity. 

The host galaxies of obscured quasars are easier to ob- 
serve and characterise since the ultraviolet and optical emis- 
sion from the central engines is obscured along the line of 
sight by intervening gas and dust. Some obscured quasars 
have been proposed to be at an evolutio nary phase where the 
host galaxies contain more gas and dust jSanders et al.l 19881 : 
lFabianiri999l ). Therefore obtaining stellar, gas and dust mass 
estimates, inferring SFR and the properties of the quasar, 
can reveal the connections between the quasar and its host, 
and test the evolution of the ratio between the black hole 
mass (M.) and the host galaxy's bulge properties, i.e. stellar 
velocity dispersion (o"*), luminosity and mass {Mbuige)- 

To derive physical properties of the molecular gas, the 
most dominant tracer of H2, carbon monoxide (CO), is ob- 
served. Observing multiple CO rotational transition lines 
provides the CO spectral energy distribution (SED), or "CO 
ladder". This can be used to infer the physical proper- 
ties, i.e. the kinetic temperature and density, of the molec- 
ular gas via the fitting of detai led models such as large 
velocity gradient (LVG) models llScoviUe fc Solomonlll974 
iGoldreich fc Kwanlll974l ). 

Other emission lines such as HCN, [CI] and [CII], 
have been observed right out to the highest redshift 
galaxies often with the help o f gravitational tensing (e.g 
ISolomon fc Vanden BoutI l2005l : IWaher et al] I2OIID . The 
atomic carbon molecule [CI] is thought to map the CO emis- 
sion, with its critical density very close to that of CO(l-O) 
10^ cm.~ ^ ). Als o, studies of the [CI] and CO in Orion by 
llkeda et al.l (|2002l ) show the emission comes from the same 
regions. [CI] can be described by the two optically-thin lines 
of a 3-level system and we can use the line ratio to directly 
determine the properties of the gas. 

The unlensed, obscured quasar, AMS12 (R.A.(J2000) 
= 17 18 22.65, dec(J2000) = -h59 01 54.3), at red- 
shift 2 = 2.767, has been observed at m ultiple 
wavelengths (iMartmez-Sansigre et all |2005| . l2006al lbl. l2009l : 
iKlockner et al.l 120091 ')." The redshift was determined from 



the Lyg, Civ and Hen lines in the optical spectrum 
l|Martmez-Sansigre et al.|[2006al '). Its mid-IR SED suggests 
the galaxy corr esponds to a progenitor of the prosont-day 
- 2L* galaxies (|Martinez-Sarisigre et al.ll2005l . l2006al lbf). A 
strong detection of the CO (3-2 ) line was first presented in 
iMartmez-Sansigre et al.l l|2009f ). prompting further investi- 
gation into the CO ladder. 

In this paper we report on the detections of two higher 
rotational transitions of CO in AMS12, along with the de- 
tections of the [CI](^Pi-^Po) and [CI](^P2 -^Pi) lines. Sec- 
tion 2 details the observations. Section 3, is dedicated to the 
dust continuum and the fit to the far-infrared (FIR) SED of 
AMS12. Section 4 details further analysis into the molecular 
gas properties using LVG models. Atomic carbon in AMS12 
is investigated further in Section 5, the results are discussed 
in Section 6. A summary of our findings is presented in Sec- 



tion 7. We assume a cosmology with Ho — 70 kms ^ Mpc ^ , 
Q.A = 0.70 and Q.m = 0.30 throughout this paper. 



2 OBSERVATIONS 

2.1 IRAM PdBI observations 

The CO(3-2) {u^cst = 345.796 GHz) transition, was ob- 
served in 2009 using IRAM's Plateau de Bure Interferom- 
eter (PdBI), in the 3 mm b and centred on 91.796 GHz 
(jMartfnez-Sansigre et al^l2009^ . The CO(5-4) and CO(7-6) 
rotational transitions (frost = 576.268 GHz and 806.652 GHz 
respectively) were observed with the PdBI 2 mm and 1.3 
mm bands centred on 152.986 and 214.148 GHz for 4 nights 
during April and May 2010. 

These observations used 6 antennas in the compact D 
configuration, with dual polarization utilising the narrow- 
band correlator. The WideX correlator was used for the 
CO(5-4) and CO(7-6) observations and covers a bandwidth 
of 3.6 GHz with a fixed spectral resolution of 2 MHz. The 
flux and phase calibrators for all observations were MWC349 
and 1637-1-574 respectively. The bandpass calibrator for the 

3 mm, 2 mm, and 1.3 mm bands was 3C345, and 3C273 
was used as the bandpass calibrator for one night of the 1.3 
mm observations. The observations were made with gen- 
erally good weather conditions (precipitable water vapour 
(PWV) ranging from 2 -8 mm), however, the 1.3 mm band 
was observed over three nights due to bad weather on those 
nights (generally conditions worsened at the end of the ob- 
serving runs with the PWV reaching > 10 mm on some 
days). 

The data wore reduced and analysed using the GILDAS 
softwar43 ( Guillotca u fc Lucas 20QO). The final data cube 
for the 3 mm band observations achieved an rms noise 
of 0.70 mjy beam~^ per 30 kms~^ channel. The result- 



ing spectrum can be seen in Figure 1(a) (as published in 
IMartmez-Sansigre et al.ll2009l ). For the 2 mm band an rms 
noise of 0.95 mJy beam~^ per 30 kms~^ bin for the final data 



cube was reached (see Figure 1(b) I. The 1.3 mm band's final 
data cube has an rms noise of 2.0 mJy beam~^ per 30 kms~^ 



bin (Figure 1(c) 



The 1.3 mm band observations revealed the [CI](^P2 — ^ 
Pi) line. This prompted an investigation into the detection 
of the [CI](^Pi Po) line in AMS12. These observations 
were carried out with the PdBI using 5 antennas in the D 
configuration over 7 days during July and August 2011. The 
observations were centred on 130.660 GHz (i^rcst = 492.161 
GHz) which falls in the 2 mm band. A total of 12.3 hours 
of on-source integration time was obtained with varying 
weather conditions (PWV ranging from 5-11 mm). Data 
were flagged on the 19 July, 02 and 05 August due to receiver 
system temperatures exceeding ~ 400 K. The uv table was 
created from the WideX correlator data with an rms noise 



of 0.577 mJy beam"^ per 20 MHz bin (see Figure |l(d) I 



The visibilities were imaged using natural weighting to 
maximise the point source sensitivity. We used the software's 
recommendations for the image and map sizes for the data- 
cube set up. The dirty images were cleaned using the Hog- 
bom algorithm down to ~ 3 times the thermal noise. The 

1 http: //www.iram.fr/IRAMFR/GILDAS 
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Figure 2. The average CO profile from the CO(3-2), CO(5-4) and 
CO(7-6) transitions. The transitions are weighted by the inverse 
square of the noises. The profile is fit with a Gaussian centred on 
-15 km s-^and with a FWHM of 290 ± 20 km s^^. 



line spectra were extracted from a single pixel centred on the 
source's position. With the resolution offered by the compact 
configuration, the source is unresolved and so extracting the 
spectra from single pixels retains most, if not all, of the in- 
formation from the source. 





Flux (mjy) 


75 


15.0 ±5 


158 


74.0 ± 20 


250 


65.0 ±6 


350 


49.4 ±5 


500 


42.5 ± 12 


1200 


3.7 ±0.6 


1400 


2.16 ± 0.19 


1960 


1.01 ± 0.06 


3260 


0.75 ± 0.09 



Table 2. The flux measurements of the FIR infrared mea- 
surements from observations by Spitzer/MIPS, Herschel/Spire, 
MAMBO and PdBI. 



2.3.1 Existing data: MAMBO and Spitzer 

The Max-Planck Millimetre Bolometer Array (MAMBO) 
observations were done at a wavelength of 1.2 mm. The 
source was observed with other sources in the AMS sample 
in blocks of typically 20 minutes. The data were reduced us- 
ing the MOPSIC p ipeline. The rms noise achieved was ~ 0.55 
mJy beam~^ (see lMartfnez-Sansigre et al.ll2009l . for details). 

We use archival Spitzer measurements of AMS 12 at 160 
and 70 fim. These were made as part of the Spitzer extra- 
galactic First Look Survey (FLS). See iFraver et al.l (120061 ) 
for details on reduction and handling. 



2.2 Results 

The linewidths for the three CO and two [CI] lines are 
~ 200 - 300 km s'V The CO(3-2), CO(5-4) and C0(7- 
6) line emission towards AMS12 are detected at IIct, 17a 
and 11a significance respectively. Figure [2] shows the stacked 
CO transition profiles with their baselines subtracted and 
weighted by the inverse squared noise. The averaged CO 
lines are well fit with a Gaussian of FWHM of 290 ± 20 
km s~^ and offset from z = 2.7668 by —15 ± 9 km s~^ 
agreeing with zco of 2.7672 ± 0.0003 as determined by the 
individual lines. 

The central frequency, flux density and full-width half 
meiximum of each line were obtained by fitting a Gaussian to 
the data, using least-squares minimisation which addition- 
ally allows an estimation of the uncertainty for the param- 
eters. These are shown in Table [T] along with other derived 
parameters. The line luminosities of both the CO transitions 
and the [CI] lines are given in both (Leo), an d K km s~^ 
pc^ (Ljo). l|Solomon. Downes fc Radfordlll992l l. 

An initial interpretation of the CO lines detected in 
AMS12, reveal they are subthermally excited, i.e. the line 
ratios are < 1. The LQQjg_4j/LQQj3_2) ratio is 0.59 ±0.05, 
and the L(Eo(7-6)/Lco{3-2) ratio 0.46 ±0.06. 



2.3 Continuum measurements 

We use a compilation of data between 70 fim and 3.0 mm to 
infer the best-fitting parameters for the FIR SED of AMS12. 
The fluxes are given in Table [2] 



2.3.2 New data: Herschel and PdBI 

The FLS field was observed as part of the Hers chel Multi- 
tiered Extragalactic Survey (HerMESfl (|01iver e t al. 201Q). 
Fluxes of the AMS objects at 250, 350 and 500 nm were 
measured off the level 2 SPIRE mosaics distributed by the 
Herschel Science Archive. Aperture photometry was car- 
ried out in apertures of 13, 17 and 23" at 250, 350 and 
500 fim, respectively, with background annuli 23-60", 30- 
100" and 40-140". Aperture corrections were derived from 
the beams associated with version 1.0 of the Spectral and 
Photometric Imaging Receiver (SPIRE) beam release note 
ISibthorpe et al.ll2010t ). sampled with 1" pixels, and the cor- 
responding beam areas applied. The beam areas assumed 
were 426, 771 and 1626 square arcseconds at 250, 350 and 
500 ^m, respectively, and aperture corrections were 1.51, 
1.52 and 1.60, respectively. 

The continuum towards AMS 12 in the PdBI 3 mm 
band was found over the line free bandwidth 684.4 MHz 
or 2053.8 kms~^, leading to an rms noise of 0.09 mJy. For 
the 2 mm band, 3.35 GHz of the line free spectrum was fit 
yielding an rms noise of 0.06 mJy. The 1.3 mm band yielded 
the continuum measurement over a spectral line free region 
spanning 2.572 GHz, giving an rms noise of 0.19 mJy. 



3 MODELLING THE FIR EMISSION 

At FIR wavelengths dust is not optically thick, and the SED 
of the radiation can be described by a modified black body. 

^ http://hermes.sussex.ac.uk 
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Figure 1. Continuum subtracted CO molecular and [CI] atomic lines from tlie PdBI 3, 2 and 1.3 mm observations, the spectra are 
in 30 km s~^ bins unless otherwise stated, (a): The 3 mm band PdBI spectrum of AMS12 with the narrow band correlator (1 GHz 
bandwidth) showing the CO(3-2) emission line, (b): Spectrum from the 2 mm band shows the CO(5-4) line, (c): Spectrum from 1.3 
mm band shows [CI](^P2 Pi) and CO(7-6) (line at km s'^). (d): The [CI](^Pi Po) hue detection, spectrum in 20 MHz bins. 
See Table [l] for the line parameters derived from the Gaussian fits to these lines. 
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-16 ±20 
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[CI](3Pi -3Po) 
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0.65 ±0.1 


11 ± 20 
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[CI](3P2-3Pi) 


214.837 ±0.016 


5.9 ± 0.9 


230 ± 30 


1.5 ±0.3 


36 ± 20 


1.7±0.3 


1.0 ±0.2 



Table 1. Observed CO and CI line parameters towards AMS12 derived from the individual Gaussian fits. ^Velocities are reported 
relative to a redshift of 2.7668. 



If the absorption coefficient of the dust, is assumed to 

follow a law oc u'^ , the emission will be given by: 



where A is a normalisation term given by: 



(1) 



A : 



+ 4)r(/? + 4) kTu' 



(2) 



The three variables are: the dust temperature Td, the emis- 
sivity index /3 and the FIR- luminosity Lfir. Here h and k 
are Planck's and Boltzmann's constants, respectively, while 
^ and r are the Riemann zeta function and the Gamma 
function, respectively. 
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In order to determine the parameters of the fit to the 
FIR SED, we use Bayes' theorem: 



P{{x\...Xn}\data) 



P{data\{xi...Xn})P{{x\...Xn}) 
P{data) 



(3) 



Here the {x\...Xn} incorporates the n variable param- 
eters of the model used. In this case, the graybody model 
with parameters T_d, /3 and Lfir fitted. The data are the 
observed fluxes. 

The posterior probabilities P[{xi...Xn\\data) may be 
split into each parameter of interest through marginalisa- 
tion. This is done by integrating the posterior probability 
density function (PDF) over the other parameters. Rewrit- 
ing Bayes' theorem with the parameters under consideration 
explicitly stated and Si, as the observable data we get, 

P(Tc, /3, LfirIS.) cx P{S,\Td,P, Lfir)P{Td)P{P)P{Lvib) 

(4) 

where the evidence is treated as a normalisation term. The 
likelihood P(5'^|ri3, /3, Lfir), is given by a Gaussian distri- 
bution: 



P(S,\TD,P,Lvm) oc e 



(5) 



with Sv,m{TD, 13, I/fir) being the predicted flux given Td, /? 
and I/FiR. Assuming Gaussian errors, maximising the like- 
lihood is equivalent to minimising the statistic where 



^(Tj,/i.£Fm)' 



X -- 

In order to get the posterior PDFs for a particular pa- 
rameter we can marginalise P(Td,P, L-pm\Sv) over the other 
two parameters. For example, to get P{Td\Sv) we integrate 
over /3 and Lfir. 



P{Td\Si,) 



D, P, LFiR\S„)dl3 dL^ 



(6) 



and similarly for P{l3\Si,) and P{LFm.\Sv)- 

shows the contours of dust temperature 



Figure 3 (a 



and emissivity index, marginalised over the FIR luminosity. 



P{Td, l3\Si,). Figures 3(b) 3(c) and 3(d) show the posterior 
distribution function for the dust temperature, emissivity 
index and Lfir marginalised over the other two parameters 
respectively. 



3.1 Far-infrared luminosity 

Figure g] shows the FIR SED for AMS12. The best fitting 
temperature is To = 88 ± 8 K and emissivity index of 
P = 0.6 ± 0.1. There are two points which are not well 
fit by the model; the 3 mm PdBI measurement, and the 
500 nm Herschel point. The 3 mm point could possibly suf- 
fer from contamination from the radio continuum from the 
AGN, AMS12 has a steep extended radio spectrum and a 
flatter, inverted spectrum within a compact 150 pc region. 
Extrapolating the radio continuum out to ~ 90 GHz us- 
ing the flatter c ompact spectra l inde x (a = —0.22 where 
Sradio OC (|Klockner et al.l [20091 ') . shows a significant 

amount of the measured flux could be due to radio emission 
from the AGN. 

The 500 /im measurement could be boosted by FIR 
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Figure 4. The graybody dust model fit to AMS12's FIR SED. 
The points are from observations from PdBI, MAMBO, Her- 
schel/Spire and Spitzer/MIPS. The best fitting temperature and 
spectral index are shown on the figure. 



emission lines. ISmail et al.l (|201ll ) have estimated the FIR 
emission lines may contribute ^ 20 — 40% to the broad- 
band flux, thus there could be possible contamination from 
the [CII] {frest ~ 158 pm) line which at z — 2.7672 is at 
595 fim, (t he Herschel/Spire 500 /xm band has a width of 
A/AA = 3. lGriflin et al]|2006h . 

Figure |4] illustrates the importance of obtaining points 
at both the longer and shorter wavelengths: it is necessary to 
probe wavelengths shorter than those where the peak of the 
emission appears, to constrain the location of this peak and 
hence the temperature. In order to constrain the emissivity 
index, long wavelength data are critical, in this case the 
data at 1.2 mm from MAMBO, and the additional longer 
wavelength data from PdBI. 

The best fitting Lfir determined from the single gray- 
body model fitting of the FIR SED is logio(LFiR/L0) = 
13.5±0.1. The Lfir can be used to determine the star forma- 
tion rates (SFRs) of the galaxies. Assuming that the Lfir 
is solely due to star formation, with young OB stars the 
main source of h eating, we can determine the SFR from the 
lKennicut3 (| 19981 ) conversion, 



SFR= 1.7 X 10" 



(7) 



In order to obtain the total SFR, the assumpt ion of an 
initial mass function (IMF) is required. Here, the ISalpeteJ 
(1955') IMF is assumed. 

The SFR for AMS12, given the Lfir determined from 
the best dust fitting model, is ~ 5300 M© yr"\ SFRs of this 
scale arc seen in only the most extreme starburst galaxies 
(see Chapman ct al. 2005; Tacconi et al. 2006; Coppi n~et al.l 
l2008l : ISolomon fc Vanden Boutiboosl . for a review). 

However, the typical temperature found for star form- 
ing galaxies in bot h the local Univer s e and at high red- 
shift, is < 50 K ( e.g iFarrah et ai1l2003l : iKovacs et aLll2006l : 
Elb az et"al1l201ir ). The temperature derived from the FIR 
SED of AMS12 is significantly higher than this. It is possible 
that a significant fraction of the Lfir is due to heating from 
the AGN, which would therefore mean the SFR determined 
here is an overestimate. 
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Figure 3. (a): The contours of the dust temperature versus the emissivity index for AMS12 marginalised over Lpia. The best fits to 
the temperature and emissivity index correspond to 88 K and 0.6. (b): The posterior PDF for the dust temperature, marginalised over 
/3 and Lfir , with best fit 89 ± 8 K. (c): The posterior PDF for the /3 parameter, marginalised over To and LpiR, with best fit value 
0.6 ± 0.1. (d): The posterior PDF of the hpjn, marginalised over Tjj and 0, the best fit is logio(Lpif(^/L0)= 13.5 ± 0.1. 



3.2 Dust mass 



The mass of the FIR d ust is also found from the Lfir by 
(e.g. iBeelen et aLll2006l ): 



Md = 



47V J, 



(8) 



where Urest is the rest frame frequency found by 
(1 + z)Vobs- The mass absorption coefficient = 
K(u^^^){i'rest/i'obs)'^ is givcu by a power law. The mass ab- 
sorption coefficient is the main source of uncertainty for 
the dust mass. We assume two different refe rence values 

of first, K{250GHz) 

= 0.04 m^kg"^ (i.e. lAlton et al.l 
l2004h. and secondly we tise K.( 24QaGHz^ = 2.64 m^kg— 1 
l|Dunne. Bales fc Edmundsll2 003^ . B^{T) is the Planck func- 
tion for a given temperature and frequency. 



3^ kT ^ 



(9) 



The dust mass for AMS12 using the fiducial values 
Td = 88 K and /3 = 0.6 is Md ~ 1.6 x lO'' Mq using 
«(250GHz) = 0.04 m^kg~\ and Md ~ 9.2 x 10^ M© using 



^(2400GHz) ~ "2.64: m kg~ . The difference between these 
two estimates highlights the uncertainty in obtaining reli- 
able dust masses. We can only estimate the dust mass to an 
order of magnitude. This is mostly due to the uncertainty in 
K and the different frequencies at which it is derived (note 
we have simply chosen two illustrative values here) , however 
other uncertainties are also large. For example the ~ 10% 
uncertainty in Td alone contributes a relative uncertainty of 
~ 30% in the dust mass. 



4 CARBON MONOXIDE 

The detection of multiple rotational transition lines of CO 
towards AMS12 gives the CO ladder. The excitation of CO 
to higher rotational levels depends upon the ambient tem- 
perature and the density of the gas. These lines are optically 
thick which makes it difficult to infer the physical properties 
of the gas directly. 
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4.1 LVG modelling 

In order to model the line intensities of the rotational transi- 
tions of CO, coupled equations of radiative transfer and sta- 
tistical equilibrium must be solved. A simplification which 
localises the problem, comes in the form of the LVG approxi- 
mation. This assumes that there are large velocity gradients 
across the area of the gas which are significantly greater 
than local thermal velocities. These act to locally trap the 
photons emitted from the de-excitation of the CO molecule 
from the J state to (J — 1). The probability that the pho- 
ton will escape the region is given by the escape probability 
which depends upon the optical depth of the transition. The 
photon trapping and the escape probability of a transition 
act to raise the intensity of the line transition above the 
background intensity (i.e. from the cosmic microwave back- 
ground radiation). 

We refer the r e ader to IScoville &: SolomonI (|l974l 'l and 
iGoldreich fc KwanI (|l974l ) for a complete derivation of the 
model in a collapsing spherical molecular cloud. 

We have used a LVG code developed and kindly pro- 
vided by C. Henkel. The LVG calculations require the col- 
lision coefficients of the molecules under consideration, in 
addition to a few input parameters. These are, an ortho-to- 
para ratio for H2, the temperature of the background CMB 
radiation, a chemical abundance of the molecule relative to 
H2 and the velocity gradient. The free parameters are the gas 
kinetic temperature, Tq, and the overall density of molecu- 
lar hydrogen, n(H2). 

We have used single-c o mpon ent LVG models using the 
collision rates from lFlowen (|200ll ). to investigate the CO ex- 
citation. In all calculations we used the H2 ortho-to-para 
ratio of 3:1, and a cosmic microwave background temper- 
ature of TcMB = 10.28 K (corresponding to the redshift 
z = 2.767). We adopted the fixed CO abundance per ve- 
locity gradient val ue of [CO]/(dv/dr) = 1 x 10~^ pc/k ms~^ 
(e.g. Wcifi, Walt er fc Scovillell2005bl : IWeifi et al.ll2007l ). 

The LVG model provides the brightness temperatures 
Tb, of each rotational transition from J = ltoJ = ll, 
which can be compared to the observed flux densities by, 



^co = ^ 



Tb 2kul 
(1 + z) c2 



(10) 



where f2 is the source solid angle (e.g. lWeifi et al ] |2007l ). Due 
to the fact that we have not resolved the source, is kept 
as a free parameter. We use the equivalent source radius ro 
which is given by ro = Da where Da is the angular 

diameter distance. 

With LVG models, there is a noted degeneracy be- 
tween the para meters Tg and n(Il2) (e.g. Wcifi ct al. 2007; 
IAo et al"]|2008l ). The degeneracy arises from the dependency 
of the line optical depth on the level populations which in 
turn depend upon the values of Tg and n(Il2). 

In order to counteract this degeneracy and constrain 
the density and temperature further, we can use the infor- 
mation we have obtained from the continuum observations. 
We present the results from the LVG modelling below, firstly 
assuming no prior knowledge on the parameters, and then 
applying a prior on the temperature from the dust analysis. 



4.2 CO LVG model results 

4-. 2.1 Assuming no prior knowledge of the gas kinetic 
temperature 

The results from the LVG modelling using flat priors for 
all the parameters, show there is no single conclusive region 
of parameter space which points to the best fitting region 
as can be seen in Figure [S] The marginalised PDFs of the 
individual parameters are also shown in Figure [S] 



The minimum temperature "floor" in Figure 5(a) 



responds to the temperature of the background radiation 
at this redshift. The minimum density "wall" reflects the 
minimum density required to excite the higher CO lines to 



the observed levels. Figure 5(a) shows there are two regions 



which are within the la contour; the region surrounding the 
triangle and the region around the square. 

Taking just our current CO measurements, we cannot 
confidently rule out either region. However, we have an up- 
per limit of ~ 15 kpc on the extent of the CO(7-6) emission 
from the unresolved PdBI 1.3 mm observations. The sizes 
which correspond to these temperatures and densities, are 
all physically possible, with the higher-temperature/lower- 
density solutions having lower sizes of the order of 1 kpc 
and less. The low-temperature/high-density solutions have 
radii of a few kpc. Figures 5(b) 5(c) and 5(d) show the 



marginalised probabilities of the unknown parameters. The 
peaks of these individual PDFs correspond to a Tg = 12 K, 
density of n(H2) — 10"^ * cm~"^ and ro — 0.7 kpc, the combi- 
nation of these individually marginalised values are marked 



by an 'X' in Figure 5(a) They do not provide a good fit to 
the data. 

The LVG model solution which corresponds to the low- 
est value has a Tg = 12 K, a density of n(H2) = lO*^'^ 
cm""^ and the size of the emitting region ro = 6.2 kpc. This 
is marked with a triangle in Figure |5(a)[ 

The low kinetic temperature of the gas from this solu- 
tion is only slightly above the temperature of the background 
radiation at this redshift. In order to get the line intensities 
above the background radiation level, this solution has a 
high gas density. The optical depths of the lines from J = 1 
to J = 7 are S> 1. At high optical depths, the photons emit- 
ted in the de-excitation of the levels remain in the region 
longer i.e. do not escape. They are able to interact further, 
driving up the number of molecules in these excited states. 

While the opacity of the lines is high, the coUisional ex- 
citation and de-excitation processes in this model dominate 
over the radiative processes and the level populations are 
in actual thermal equilibrium. The excitation temperatures 
of these line transitions is equal to the kinetic gas tempera- 
ture since the system is in local thermodynamic equilibrium 
(LTE). 

At higher transitions (J ^ 8) the line optical depths are 
< 1 and the system is no longer in LTE. The lines undergo 
coUisional de-excitation and the intensities quickly decline to 
levels no longer detectable above the background radiation. 
This manifests itself in the sharp decline of the solid line in 
Figure [6] above J — 7. 

The dashed line in this figure is the CO SED corre- 
sponding to a secondary minimum region within the 

with Tg = 135 K, n(H2) = 



Icr contours of Figure 5(a 

^'^ "m~^ and ro = 0.7 kpc (the square in Figure 5(a) I 
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Figure 5. (a): The contours of the temperature and density PDFs marginalised over the size. The open triangle marks the model of 
least x^i while the open square marks an arbitrary spot in the higher temperature-lower density/size region, within the la confidence 
interval. The best values of the marginalised individual parameters are marked with an 'X'. (b): The PDF of the gas temperature 
marginalised over the density and size parameters, (c): The PDF of the gas density parameter marginalised over the temperature and 
the size. (d):The PDF of the size of the emitting region of gas marginalised over the temperature and density. 



Here both the radiative and collisional excitation and de- 
excitation have significant effects on the line intensities and 
the full statistical equilibrium analysis must be considered. 

Figure [6] illustrates the difference between this high- 
temperature/lower-density solution (the dashed line), and 
the previous low-temperature/high-density solution (solid 
line) is most pronounced at the higher J levels. Observations 
of the J ^ 8 transitions are needed to further constrain the 
LVG modelling at these higher transitions. 

High resolution imaging of the gas (and dust) in AMS12 
would provide a measurement on the size and would con- 
strain the model further. 



4-2.2 Using the dust temperature PDF as a prior for To 

Here we work under the assumption that the gas and dust 
arise from the same regions, and are therefore at the same 
temperature. Indeed dust shields the gas from the ultravio- 
let and optical radiation preventing it from dissociating the 
molecules. 

Gas and dust studies in other high redshift galaxies have 
shown both the dust and gas are compact on scales of less 



than ~ 4 — 8 kpc suppo rting the as sumption that they arise 
from the same region llXacconi et al. 200f^; iRiechers et al.l 
l2006l : IWeifi et al.ll2007l : [Younger et al.,.200a ). The dust tem- 
perature is often used to either constrain or ju stify a par- 
ticula r gas kine tic temp erature (s ee for example ^We ifi et alj 
,200i lAo et~alli2008 ; Greve et aL|[200a . among others). 

Using the PDF of the dust temperature as the prior for 



the gas temperature, (i.e Figure 3(b) where Td peaks at 89 
K), P(Tg) = P(Td) in Equation 31 we can rule out the lower 
temperature regions. With the prior on the temperature dis- 



tribution, the contour map in Figure 7(a) shows a tight con- 



vergence on the temperature and density. The values of the 
least in this model are approximately equivalent to the 
peaks of the marginalised PDFs of the parameters which are 
shown in Figure [T] Due to the small number of data points 
and the shapes of the likelihood and prior PDF, the prior is 
having a greater effect than the likelihood on the PDF. The 
values for the kinetic temperature, density of the gas and 
the size of the emission region are 90 ±8 K, 10'^'^='=°-^ cm"'', 
and 0.8 ± 0.04 kpc, respectively. 

The LVG model corresponding to these values of tem- 
perature, density and emitting region as determined with 
the use of the dust temperature PDF as a prior for Tq, 
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Figure 7. (a): The 1, 2, and 3o" contours of the temperature and density PDFs marginahsed over the size. The dust temperature is 
used as a prior for the gas temperature. The best values of marginahsed temperature, and density are shown in each figure marked with 
an 'X'. (b): The marginalised PDF of the gas temperature using the prior of P(Td)=P(Tg). Best fit with Tq = 90 ± 8 K. (c): The 
marginalised PDF of the gas density parameter logio[n(H2)/cm~'^]=3.9±0.1. (d): The marginalised PDF of the size of the emitting 
region of gas best fit with logio[ro/pc]=2.9±0.02. 



gives the CO ladder displayed in Figure [8] The combina- 
tion of these parameter values fall within the Icr contours of 



Figure 5(a) making them an acceptable fit to the CO SED 



Eissuming no prior knowledge. 



5 ATOMIC CARBON 

The [CI](^P2 -'^ Pi) hue was detected in AMS12 to an 7(t 
significance which prompted a search for the [CI](''Pi — ''Po) 
line. The [CIK^Pi Po) line was detected to an 8cr signifi- 
cance, these lines can be seen in Figure[T]and the line proper- 
ties are given in Table [T] This is the first unlensed, high red- 
shift galaxy with both the [CI](^Pi-^Po) and [CI](^P2-^Pi) 
lines detected. 

With both the upper and lower fine structure atomic 
carbon lines detected, we may directly determine physical 
properties of [CI] in AMS12. 



5.1 Excitation temperature and mass of atomic 
carbon 

The excitation temperature of [CI] can be directly deter- 
mined, assuming the lines are optically thin, from the rati o 
of L[^il(3p2-3pj) to L[^i](3pj_3p^,) (e.g. IStutzki et al.lll997l'). 
In order to relate the excitation temperature of [CI] to the 
gas kinetic temperature, we require the [CI] excitation to 
be in LTE (much in the same way we can estimate the Tq 
directly from CO in LTE). If this is not the case, the exci- 
tation temperature of [CI] could be lower than the kinetic 
temperature also. 

The line column densities may be found using the in- 
tegr ated brightness temper atures of the lines (see Appendix 
A of ISchneider et aLll2003l . for a complete derivation). The 
excitation temperature is given by the ratio of column den- 
sities expressed by the ratio of the statistical weights of the 
levels and the Boltzmann factor. 
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Figure 6. The CO SED with the SED given by the LVG model 
giving the lowest with Tq = 12 K, n(H2) = 10^-^ cm~^, and 
the CO region size of tq = 6.2 kpc (solid line). This corresponds to 
the open triangle in Figure [5( a) | A second line, showing a model 
from within the higher temperature-lower density/size region (Tq 
= 135 K, n(H2) = 10^'* cm~^ and ro = 0.7 kpc), corresponding 
to the open square in Figure [5(a)| is there for comparison (dashed 
line). 
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Figure 8. The CO SED with the SED given by the LVG model 
with Tg = 90 K, n(H2) = lO^'^ cm^^^ and the CO region size of 
ro = 0.8 kpc, these have the highest probabilities once the prior 
on the temperature, P(To), is applied. 



rearranging and equating constants the expression for the 
excitation temperature is found to be, 



Tex — 



38.8 



(12) 



where _Rci is the ratio of the line brightness temperatures, 
these may be replaced with the L"'^ of the hues as L""" oc Tf,, 

RCI = L[CI](3p2_3p^)/L[(.ij(3p^_3p^). 

The Tex can t hen be used t o find the [CI] total column 
density and mass. IWeifi et all (|2003l ) derive the beam av- 
eraged [CI] column density in the optically thin limit and 
use this and the area of the emitting region (given by the 
sohd angle subtended by the source convolved with the beam 



multiplied by the angu lar distance squared, Q,a*bDA) , to de - 
rive the mass of [CI]. ISolomon. Downes fc RadfordI (|l992l ) 
express the line brightness luminosity related to the emit- 
ting area, = 23.5f2s«6-D_t Jco(l + , where the lumi- 
nosity distance is related to the angular di stance via Da ~ 



DL/{l+zf. Thus. lWeifi et al.l ll2003l.l2005al). use this expres- 
sion o f and Equation 3 in ISolomon. Downes fc RadfordI 
(1991) in order to determine the mass of [CI] in an unre- 
solved source using L""" 



[CI]{3Pi-3Po) 



[Ma] ^ 

5.706 X 10"''Q(Te 



T T 



-3Po) 



K kms^i pc^ 



(13) 



where Q(Tcx) = 1 + 3e~'^i'"^'=- + Se-'^^Z'^" 
function for [CI]. Ti = 23.6 K and T2 = 
energies above the ground state. 



is the partition 
62.5 K are the 



5.2 The [CI] temperature and mass towards 
AMS12 

The relationships defined above require the [CI] lines to 
be optically thin. We can test this requirement assuming 
a Tex and size (se e equations A6 and A7 in the appendix 
of ISchneider et all [2003 ). Firstly, assuming Tex = Tq and 
a size of 0.8 kpc (i.e. from the LVG model) we infer line 
optical depths of ~ 0.3 for both the [CI](^Pi Po) and 
[CI](^P2 — -Pi) lines. If we assume a larger size (i.e. 2 kpc), 
the optical depths are ^ 0.05. 

From the line luminosities of the atomic carbon in 
AMS12 (see Table [!}, we determine the line ratio Rci ~ 
0.85 ± 0.24 and hence a [CI] excitation temperature of 
42.7 ± 10 K. This yields a [CI] mass of (1.5 ± 0.5) x 10^ 
Me. 

The [CI] excitation temperature determined does not 



fall within either la temperature regions in Figure 5(a) 



Thus when we fit the LVG model with no prior on the tem- 
perature, the Tex of [CI] is not in agreement with the gas 
temperatures indicated by high-excitation CO. It is in even 
less agreement with the gas kinetic temperature we deter- 
mine using the dust temperature as a prior. This could be ei- 
ther suggesting that the [CI] is not in LTE. Or alternatively, 
the [CI] emission arises from a more spatially-extended, 
cooler molecular gas component than the gas probed by the 
higher-excitation CO lines. 

The optical depth of the [CI] increases with decreasing 
temperature, using Tex ~ 43 K the line optical depths are 
calculated to be ~ 1 when assuming the compact size of 0.8 
kpc. However, assuming the size is instead more extended for 
example, 2 kpc, the line optical depths are ~ 0.1. It is likely 
that our assumption of optically-thin lines is appropriate 
as even with the more compact size we are at the limit of 
optically-thin lines, and the possibility that the [CI] emission 
region is more extended would lower the line optical depths. 

The LVG model analysis of the high-excitation gas has 
indicated this to be relatively compact with a radius ~ 1 
kpc. The low-excitation temperature of [CI] which in other 
high redshift sources is broadly in agreement with the dust 
temperature (|Walter et al.|[201ll ). could be alluding to a sec- 
ond, cooler, more extended region of gas. Since we expect the 



Gas and dust in AMS12 11 



CO (1-0) to trace the same region as the [CI] emission, if we 
were to observe CO (1-0) we could test if this low-excitation 
gas component exists in AMS12. 

This has been seen in high redshift submillimetre galax- 
ies (SMGs). High resolution imaging of CO(l-O) has re- 
vealed the CO(l-O) to be more e xtended than the higher 
J emission ( i e. CO(3-2) or (4-3)) jlvison et al.ll2010l . I2OIII : 
ICarilli et all I2OI0I : iRiechers etall l2011bt ). While resolved 
CO(l-O) in strongly lensed quasar hosts shows the CO(l-O) 
to be compact and similar to the CO(3-2) emission, assum- 
ing the low-excitation gas has th e same ma gnification factor 
as the high-excitation gas (Ric chers et al.L 2006. 2011a). 

[Cerin k, Phillips! (|2000l ) observed a relationship be- 
tween Lqo(i-o) Ljci](i-o) from a survey of low redshift 
galaxies; Lj;^i](i_o) = 0.2±0.2 L^r,(i pj. lWalter et all (|201ll ') 
observed that their sample of (mostly l ensed), high red- 
shift s ources support the relationship from lCerin fc Phillips! 
l|2000l ). For the majority of their sources, CO(3-2) was the 
lowest transition observed, in order to get an estimate on 
CO(l-O) they used Lco(i_o) ~ 0.9 Lco(3-2) ■ From this, they 
determined that LJ^jj^j^.qj = 0.29±0.13 Lco(i-o)- 

However, in the case of one SMC (SMM 
J163658+4105), which llvison et all (120111 ) observed in 
CO (1-0), the ratio of the brightness temperature lumi- 
nosities of CO (3-2) and CO (1-0) is actually 0.54 ± 0.12 
signifi cantly lower than the assumed 0.9 that IWalter et al 



2011) used. Using the CO(l-O) strength from llvison et al 
20111 ). the ratio of Lco(i-o)/Ljci](i-o) then becomes 0.14, 
in disagreement w ith the high redshift relationship from 
I Walter et all l|201ll ) and closer to the low redshift relation- 
ship. The ratio of Lco(i-o)/L[ci](i-o) high redshift might 
well be similar to that at low redshift, once differential 
magnification of high- and low- J CO lines in lensed sources, 
has been taken into account (e.g. Deane et al. 2012b, in 
prep.). 



these methods presented here assume that the molecular gas 
is represented by a single highly-excited component. We can 
also use the [CI] ( ^Pi — Po) line and the relationship from 
ICerin fc Phillips! (|2000h to estimate the CO (1-0) brightness 
temperature luminosity without this assumption. 



6.1.1 Method A 

Assuming the transitions lower than CO (3-2) are ther- 
malised, we can also derive a total gas mass, M(H2), from 
using the relationship in Equation 1141 
Under the assumption that CO(3-2) is thermalised, 
Lco(3-2) = Lco(i-o) = (4.2 ± 0.4) x 10^° K km s'^ pc^ 
This gives a gas mass of M(H2) ~ 3.3 X loi°((f^)M0. 



-■^00(1-0) 



6.1.2 Method B 

Method B uses the best fit LVG model shown in Figure [6] 
(solid line. To = 12 K, n(H2) = lO'^'^ cm"'*). The CO(l-O) 
line brightness temperature Tb is given by the LVG model, 
and following the use of Equation [10] to convert to observed 
flux, we calculate Lco(i-o) ~ 4. 8 x lO'^" K km s~^ pc^. 
The gas mass determined from Method B is M(H2) ^ 3.9 x 
lOi°(^)M0. 



6.1.3 Method C 

We have used the dust information we have for this source 
and applied the prior distribution for the gas temperature 
P(Tg) = P(Td). The best fit model is shown in Figure [8] 
(Tg = 90 K, n(H2) = 10^'^ cm"^). We can again use this 
LVG model's CO (1-0) Tb and convert it to a luminosity of 
Lco{i-o) ~ 3.2 X 10^" K km s"^ pc^ This gives Mu^ ~ 
2.6 X lO"'(ji|)M0. 



6 DISCUSSION 
6.1 Gas mass 

From the measured brightness temperature luminosity of the 
CO gas (LSo),the H2 mass can be found using the relation: 



j^ M(H2) 



Mg 



Mq{K kms-i pc2)- 



T T 



(14) 



(see ISolomon fc Vanden Boutll2005l). where M(H 2 ) incl udes 



He and is therefore « Mgas. Downes fc SolomonI (|l998h de- 



termined the constant a empirically from a study of a sam- 
ple of ultra-luminous infrared galaxies (ULlRGs) and high 
redshift galaxies to be a = 0.8. We assume this value 
of a = 0.8 hereafter, although we note that the ac- 
ceptable values are in the range 0.3 - 1.3 (see Table 9 of 
iDownes fc Solomon 1998). 

The brightness temperature luminosity of CO is best 
represented by Lqoji.q)- This requires a measurement of 
CO(l-O) that we do not have, but we can calculate esti- 
mates of the CO (1-0) fiux via various methods. Three of 



6.1.i Method D 

Using the relationship relating ^^coi^-o) to L[^jj(^_g) from 
.Ger in fc Phillips! (|2000! 'l. we estimate the Lco(i-o) = 6.1 x 
10^" K km s"^ pc^ We use the low redshift relation- 
ship as the results from !Walter et all (|201l! ) are estimated 
from the higher CO (3-2) transition and not CO (1-0). The 
gas mass inferred from this CO (1-0) strength is Mu^ ~ 
4.9 X lO^''(j^)M0. This is larger than the masses determined 
from the previous three methods which is expected if the [CI] 
observations are revealing a low-excitation, more diffuse re- 
gion of gas, i.e. such as the low-excitation gas re servoirs seen 
in the CO (1-0) observations of some SMGs (llvison et al.l 
l2010l . l201ll : ICarilh et al.ll2010l : IRiechers et al.ll2011bl ). 

The results are displayed in Table [H] The first three 
methods {A to C), assume that the high- J CO lines trace 
the total molecular gas. The variance between the mass es- 
timates from these three methods yields an estimate of the 
uncertainty in the mean value (Mhs) = (3.3±0.9) x 10^" {f^) 
M0. However, with both [CI] lines detected we have found 
a [CI] excitation temperature which does not agree with ei- 
ther of the LVG model solutions. It is thought that [CI] 
traces the low- excitation CO, as the critical densities of [CI] 
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and CO (1-0) are very similar and there is observational ev- 
idence of [CI] an d CO arising from the same regions (e.g. 
Ilkeda et al.ll2002l ). The assumption that the high-excitation 
CO lines trace CO (1-0) may be incorrect and lead us to 
underestimate the total gas mass. 

Indeed, studies of high redshift SMGs, show the resolved 
CO(l-O) evident ly arise frorn exte nded^ low-excitatio n gas 
(see for example llvison et al.ll20ld. 12 011: Caril h et al.ll20ld . 
and references therein! . Ilvison et al.l "( 2011 ) found that using 
CO (1-0) to determine the gas mass in four SMGs at z ^ 2 
gave masses ~ 2 times higher than masses determined from 
the CO(3-2) or higher J lines, i.e. the assumption that the 
higher J transition (CO (3-2) or CO (4-3) in most cases), is 
thermalised was incorrect. 

Single component LVG models of SMG s have been 
show n to underestimate the CO(l-O) line fi.e. ICarilli et al.) 
I2OI0I : iRiechers et al. 2011b. ). Direct observations of the 
CO(l-O) line in AMS12 are needed to test whether a sig- 
nificant component of low-excitation gas is present. 

Hence, from here on we shall use the gas mass de- 
termined from the CO (1-0) strength estimated from the 
[CI](^Pi -^Po) line (method D), Mh^ ~ 4.9 x 1O"'(^)M0, 
although we acknowledge there is considerable uncer- 
tainty in this estimate due to the large scatter in the 



sider the LVG model which gives the lowest x value, Tq 



12 K, n(H2) 



10"= 



cm and ro = 6.2 kpc, and as- 



^co (i-o)/L[ci](i-o) 



relationship. 



6.2 Discriminating between LVG models 

We discuss various methods we use to discriminate between 
the LVG models. First, we consider the models' outputs and 
the constraints future observations can make, secondly, by 
estimating the mass given by the models assuming the gas 
is distributed in a thin disc, we compare these with the gas 
mass from Section FB.ll Then we explore the possibility of a 
two component LVG model inspired by the excitation tem- 
perature of [CI]. 



6.2.1 From the LVG models' output line brightness 
temperatures 

The sharp decline in the intensity of the higher- J tran- 
sitions in the low-temperature/high-density model can be 
used with future observations to discriminate between the 
models. Using the output Tb's we can predict the possible 
line strengths of the higher J transitions, particularly us- 
ing the CO (8-7) and CO (9-8) transitions which happen to 
be in observable windows. These would be 500(8-7) ~ 1-1 
mjy and Soo(9-8) = 0.03 mjy respectively (see Figure [6]). 
However, if we instead use the Tb's from the solution where 
we used the dust temperature as a prior for the gas tem- 
perature the line strengths are Sco(a-7) ~ 5.9 mJy and 
500(9-8) =3.1 mJy respectively (see Figure [S]). 

The significant differences in these line strengths mean 
that observations of higher CO transitions in this object, 
would conclusively constrain the region of parameter space 
of the LVG model which describes the conditions of the gas. 



6.2.2 From mass estimates using volume arguments 

Another way we may be able to rule out one of the models 
is by considering the spatial distribution of the gas. Con- 



sume the molecular gas is distributed in a thin disc with 
height H, of uniform density. We can calculate the vol- 
ume of the gas given the density. We use the findings of 
iDownes fc SolomonI (|l998l l that the molecular gas in the cen- 
tral regions of ULIRGs is not a collection of separate clouds 
undergoing self-gravitation, but rather clouds fused together 
to form a disc of more or less constant density. They mod- 
elled the structure of these gas discs and found the average 
height of the discs to be _H" ~ 58 pc. Taking this value as 
the height of the assumed disc of gas for AMS12, we can 
calculate the volume of the disc given ro = 6.2 kpc. 

From this volume, and the best fitting density of the 
gas n(H2) — lO"'^ cm~'^, we can estimate the gas mass in 
this volume to be ~ 4.1 x 10^* Mq. This is four orders 
of magnitude higher than the gas mass derived from our 
observations (see Section [6.1|) . Table O shows the gas mass 
Mh2, for the four methods described in Section r6.ll 

On the other hand, the LVG model giving the parame- 
ters To = 90 K, n(H2) = 10^ ® cm"^ and ro = 0.8 kpc, gives 
a mass of 4.6 x lO^'' M© assuming the same disc height of 58 
pc. This is very similar to the masses in Table [3] determined 
by Equation [TH 

Though we note we have made various assumptions to 
estimate the gas masses in Section [6.11 this crude mass es- 
timate from the volume argument, distinguishes the LVG 
model solutions from each other. The estimate from the vol- 
ume of the low-temperature/high-density solution does not 
agree with the gas masses determined in Section [6.1l while 
the mass from the solution using the dust temperature as a 
prior, does agree. 



6.2.3 Two-component model 

IRiechers et"all l|2011bl ) found in the two SMGs they studied, 
the CO ladders were best fit with two component LVG mod- 
els; a dense, high-excitation component to fit the higher-J 
CO lines, and another low-excitation component to fit the 
low-J CO measurements. 

To investigate the possibility the [CI] gas in AMS12 is 
from a lower-temperature, more diffuse region, we fit a sec- 
ond component to the LVG model. We used the CO (1-0) 
estimate from the [CI]('^Pi — Po) line, and for the high- 
excitation component we keep the previous LVG model so- 
lution using the dust prior (from Figure [8|. 

We fix the temperature of the second component at 
42.7 K (the [CI] excitation temperature), and tie both com- 
ponents to the estimate of the CO (1-0) strength. Since we 
only have one data point, we cannot fit for both the size and 
the density at the same time, so we fit each independently 
keeping the other fixed at a fiducial value. Firstly, we keep 
the density fixed to « 10^'^ cm~^ and fit for the size. This 
density is chosen as a representative of a diffuse component 
and corresponds to the critical density of CO (1-0). The best 
fit to the radius is 2.8 kpc and the resulting LVG model is 



shown in Figure 9(a] 

Next, we fix the size of the low- excitation gas, while 
continuing to hold the temperature at 42.7 K, and fit for 
the density. We choose two arbitrary sizes; 1 kpc, which is 
similar to the best fitting size of the single component LVG 
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Method 


CO(l-O)/ " 
[K kms-i pc2] 


MhJIQIO 

[Mq] 


'^[CI1(1-0)/-'^CO(1-0) 


A 


4.181 


3.34 


0.29 


B 


4.837 


3.87 


0.25 


C 


3.246 


2.60 


0.37 


D 


6.085 


4.89 


0.20* 



Table 3. The calculated CO(l-O) line strength determined four ways. A: from the assumption that CO(3-2) is thermalised; 
B: from the unconstrained best fitt ing LVG model; C: from the LVG model assuming P(Td) = P(Tg); and D: from the 
[CI](3Pi -3 Po) strength using the ICerin fc Phillips! ||200C|) rela tionship. The [CI]/CO relationships for each CO(l-O) are 
determined and displayed. * This is the lGerhT^ Phillipj l(2000l 'l relationship that we have assumed to derive the CO (1-0) 
mass. 



model solution to the high-J CO lines, and 6 kpc, which is 
significantly more extended. For the second component so- 
lution with a fixed size of 1 kpc, the best fitting density is 



n(H2 



10 cm . The two component LVG model (com- 



prising of this solution and the best fitting solution of the 
high-J CO lines with the dust temperature prior applied), 
is shown as the solid line in Figure 9(b) It does not provide 
a good fit to the observations. 

With the second component size fixed at 6 kpc, (and 
temperature remaining at 42.7 K), the best fitting density 
is n(H2) = 10^'^ cm~^. The corresponding two component 



LVG model is shown in Figure 9(c) This second component 
solution provides a better fit to the observed CO lines. The 
estimates shown in Figures [9(i[)] and[9(^ are consistent with 
a more diffuse and extended gas which is only detectable via 
the low-J CO lines and the [CI] lines, while still providing 
good fits to the high-J observations. 

We note that by construction, our method will only 
select solutions with weak emission of the high-J lines. Since 
we have begun by fitting CO(3-2), CO(5-4) and CO(7-6) 
with a single component model, the extra component used 
to fit CO(l-O) must produce negligible fiux in the CO(3-2), 
00(5-4) and CO(7-6) lines, otherwise their predicted fiuxes 
will be higher than those observed. 



6.3 Dynamical mass 

We can estimate the dynamical mass of the system from the 
FWHM of the CO lines, assuming a characteristic radius via 
iNeri et al.ll20o3 ): 



MdynSin'^iX 4, ( AVfwhm 

' = 4 X 10 ' 



\ kms~i 



Studies have shown the virial mass estimate is reason- 
able even if the gas is clumpy (see [Pa ddi ct al. 2010). 

We can use the CO linewidths to estimate the dynami- 
cal mass, however, if the [CI] traces a low-excitation region 
of gas, we should use the [CI] linewidths. Given the large 
uncertainties in the linewidths in Table [T] the CO and [CI] 
linewidths are very similar and agree within 2a. Using both 
the [CI] linewidths, we use AVfwhm ~ 260 km s~\ This 
will give an approximate estimate on the dynamical mass. 
Assuming two arbitrary radii encompassing a range of sizes, 
i.e. the radii from Section [6.2.31 r\ = 1 kpc and r2 = 6 kpc. 



we can estimate two values of the dynamical mass of AMS12. 
Using ri, the dynamical mass is AfdynSin^i ~ 2.7 x 10^ M©, 
while using vz, gives MdynSin i ^ 1.6 X 10^'' Mq. 

The dynamics are dominated by the molecular gas 
and stars, with the dark mat ter and ionized hydrogen sub- 
dominant (|Daddi et all 12010). The steUar mass in AMS12 
is fou nd from the mid-/ncar- infrared SED by iLacv et al.l 
(|201ll ). to be 3 X 10" Mq. The dynamical mass es- 

timates make it difficult to account for the stellar mass. It 
may be that the radius estimates we have used are not ad- 
equate to encompass all the stellar mass. The constraints 
on the inclination angle assuming the radii above are severe 
when considering the stellar and gas mass, for example, if 
ri = 1 kpc, then the inclination angle 5°. While using 
the radius estimate r2 = 6 kpc, the constraint on the incli- 
nation angle relaxes slightly to i ^ 13°. These arguments 
suggest the host galaxy is seen face-on. 

The radio spectrum o f AMS1 2 has been inve s tigate d 
in iMartfnez-Sansigre et all (|2006bl ): iKTockner et al.l (|2009l ). 
This object has a steep extended radio spectrum, and nar- 
row optical emission lines pointing to torus obscuration, i.e. 
the orientation of the central engine and its obscuring mate- 
rial is closer to edge-on to the observer. The low inclination 
angle of the host from the dynamical mass estimates, sug- 
gest the central AGN region and the galaxy's stellar, gas 
and dust regions are not aligned. Together with the ob- 
served narrow emission lines, this suggests AMS12 is ob- 
scured by the torus and not by dust in the host galaxy (see 



scurea by tne torus and not by i 
iMartfnez-Sansigre et aLll2006al ). 



The caveats to using this dynamical mass estimate are 
many. The high-excitation CO lines may be tracing a sep- 
arate gas component than the [CI] lines, and therefore the 
radius estimates based on the high-excitation CO LVG mod- 
elling are tenuous for the possible low- excitation component. 
Though the [CI] linewidths are similar to the CO linewidths, 
this does not immediately place the [CI] at the same region 
of the high-excitation CO, i.e. if [CI] is alluding to a more 
massive while more extended low- excitation gas reservoir. 
Resolving the detected CO and [CI] lines would constrain 
the sizes of these emitting regions, which may be possible 
with the upcoming IRAM Northern Extended Millimeter 
Array (NoEMA) upgrade to PdBI. In particular though, re- 
solved CO(l-O) in this object would significantly improve 
the dynamical mass estimate in this object. 



14 Schumacher et al. 



10r-r 




2 - 



n l , , , I I , '~i I J 

1 23456789 10 11 
Rotational Quantum Number J„pp„ 

(b) 



10|— r 




1 23456789 10 11 
Rotational Quantum Number J„pp„ 

(c) 

Figure 9. Two component LVG models fitted to the CO obser- 
vations and an estimated CO(l-O) line strength of 0.57 mjy. The 
low-excitation component is given by the dashed line, while the 
high-excitation component (from Figure [Sjl is given by the dot- 
dot-dash line. The combined model is given by the solid line, (a): 
The low-excitation component with Tq fa 43 K, n(H2) = 10^'^ 
cm~'^ and size of rg = 2.8 kpc. (b): The low-excitation compo- 
nent in this figure has the solution Tq = 42.7 K, n(H2) = lO'^'^ 
cm""^ and a fixed size of ro = 1 kpc. (c): The low-excitation com- 
ponent here has Tq 43 K, n(H2) = 10^'^ cm~'^ and size of rg 
= 6 kpc. 



6.4 [CI] abundance and cooling contribution in 
AMS12 

The abundance of [CI] in AMS12 can be determined by 
X[CI]/X[H2] = Mci/(6MhJ, where Myci] ~ 1.5 x 10^ Mq. 
The abundance assuming the high-excitation gas can be 
used to estimate the gas mass is 7.8 x 10~^. While if there is 
a lower temperature, more diffuse component of the gas, 
the would be higher. The Mhj estimated from the 

[CI](3Pi-3Po) line is « 5x10^° Mq , giving a [CI] abundance 
of X[CI]/X[H2] = 5.2 X 10-^ This indicates the molecular 
gas i s already enri c hed a t this redshift, supporting findings 
from I Walter et af] (|201ll ) . 

The ratio of L[ci]{i_o)/Lfir provides a measure of the 
cooling contribution of [CI]. For AMS12, L[ci]{i-o)/Lfir = 
1.5><102lL.wliich appe ars to be typical for the quasar sources 
of I Walter et"aLl (|201ll ). Their sources are split into quasars 
and SMGs, and while quasars have L[ci]{i-o)/Lfih, ratios 
similar to that of AMS12, the SMGs ratios are around an 
order of magnitude higher. This could be due to the AGN 
contribution to the Lfir in the quasars, as discussed earlier. 
Overall, the [CI] lines are not major coolants, in fact they are 
negligible compared to the cooling by the dust continuum. 

6.5 AGN bolometric luminosity and the scale of 
AGN heating 

From our current observations of the dust and high- 
excitation gas in AMS12 the temperatures we determine 
are shown to be beyond the capabilities of heating by star 
formation alone. Empirical observational evidence and ra- 
diative models of star-forming galaxies show that the dust 
temperatures reach ^ 50 K (see fo r example iKovacs et al.l 
I2OO6I : ISiebenmorgen fc Kriigelll2007l) . The typical dust tem- 
peratures derived from FIR SED fitting in studies of high 
redshift star-forming gala xies, range from 30-60 K with the 
average being ~ 35 K (e .g. iKovacs et allbood : [ Coppin et al.l 
I2OO8I : lElbaz et al.ll2011^ . These temperatures are typical of 
local galaxies where heating of the d ust is dominat ed by 
young stars (e.g. Farrah ot al. 2003; E lbaz et aLlbOllI ). 

However, in the hosts of luminous AGN, the FIR emis- 
sion could also be heated by the AGN. APM 08729-1-5255, 
F10214, BR 1202-0725 and Cloverleaf for example, have hot- 
ter FIR dust than typical star-forming galaxies and other 
SMGs. Detailed studies of the dust and gas in these ob- 
jects have revealed the dust is compac t which supports the 
possibility of significant AGN heating (ISolomon et al.|[2003l : 
iRiechers et al1l2006l : IWeifi et al.|[2007l : lAo et al.ll20oi ). 

We can estimate the scale of the heating from the 
AGN using the bolometric luminosity of AMS12 from the 
broad-band data between 3 .6 and 24 /im (as done by 
iMartmez-Sansigre et al1l2009l '). For AMS12, an AGN of Lboi 
= 2 X 10^'^ Lq, we assume that Lt,„ is ^ 0.25I/boi, and using 
the FIR SED fitted parameters {To = 88 K and /3 = 0.6), 
we find the scale of AGN heated dust out to 88 K is 2.7 kpc. 
This result is obtained assuming that the ultraviolet pho- 
tons travel unhindered up until the radius where the dust 
becomes self-shielding (i.e. at large optical depths), thus 2.7 
kpc may be thought of as the characteristic radius to which 
the dust is he ated to this temperature by the AGN (see 
lBarvaini3ll987l . for details). Clearly, since the characteristic 
scale of the dust is indeed found to be 2 kpc in many ob- 
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iects jGreve et al.ll2005l : iTacconi et aLll2006l : I Younger et"all 
l2008l 'l. the dust temperature observed in AMS12, Td = 88 
K, could be achieved through heating from the AGN. 

Note, this estimate of the bolometric luminosity does 
not take the -Lfir into account, and therefore may be con- 
sidered a lower limit. If we account for the Lfir which we 
believe is also attributed to the AGN (3 x 10^^ L©), our 
new estimate of the Lboi is 5 x 10" Lq for AMS12, which 
could heat the dust to 88 K out to around ~ 4 kpc. 



6.6 Black hole, stellar, gas and dust masses 

We estimate the dust mass to be between 9.2 x 10^ M© and 
1.6 X 10^ M0 using the Lfir and two different mass ab- 
sorption coefficients (this range is typical for high redshift 
sources, see for example ISolomon fc Vanden BoutI (|2005h '). 
The dust masses vary greatly due to the process of extrap- 
olating the mass absorption coefficients to the rest frequen- 
cies, a power law which depends upon /3. The gas mass from 
the [CI] observations is determined to be 4.9 x 10^" M©. 
This yields a gas-to-dust mass ratio between ~ 30 — 530. 

We can estimate the Eddington limited black hole mass 
of this system given the Lboi. If accreting at ^ 100% of the 
Ed dington rate (reasonable for quasars at high redshifts, 
see iMcLure fc Dunlodl20oA . the black hole mass estimated 
from the Lboi determined from the mid-infrared SED is 
M. > 6 X 10® M0. The revised Lboi (including the Lfir), 
leads to an estimate of the black hole mass M, ^ 1.5 x 10^ 
Mq . We have assumed that all of the Lfir is attributed to 
heating from the AGN, though we acknowledge that star 
formation is likely to contribute, we do not have a measure 
of the extent of this contribution. 

The stellar mass of AMS12 has been estimated by 
iLacv et al.l (|201ll ) to be M* = (3.2 ±0.3) x lO" Mq. We as- 
sume this stellar mass is located entirely within the bulge as 
AMS12 is probably the progenitor of a present-day elliptical 
galaxy. This means AMS12 has already the stellar mass of 
a present-day 2L* galaxy. 

Given the gas mass of A/h2 ~ 4.9 x 10^'' L© we derive 
from the CO observations, were this to be converted into 
stars with 100% efficiency, it would still only increase the 
stellar mass by ~ 15%. Observations of CO(l-O) are needed 
to probe the lower-excitation gas and give a more accurate 
value of the gas mass. However, it seems unlikely that the 
stellar mass will increase significantly, unless we have under- 
estimated the gas mass by a factor of ~ 10. 

Assuming M* = Altuige, we have a M, /Altuige ratio 
for AMS12 of > 0.005. This is significantly higher than 
th e relationship determine d fr om nearby (z ~ 0) g alaxies 
bv lMarconi fc Iluntl l|2003l ) and lHaring fc Rixl l|2004h where 
M./Mby,ige ~ 0.002 with a scatter of ~ 0.3 dex. 

The M, — Mfjuige relationship has been i nvestigated at 
higher redshifts up to z ~ 4. For example, iMcLure et al.l 
(|2006h investigated the rela tionship in radio-loud galaxies at 
z ~ 2. iDecarli et"ai1 ll2qiOl) stu died a sample of 96 quasars 



out to 2 ~ 3. |Peng et al.l 12006a,'?) used both gravitationally 
lensed and non-l ensed galaxies to study the r e lation ship out 
to 2 ~ 4.5, while iTargett. Dunlop fc McLurj (|2012h studied 
2 ~ 4 quasars. The high redshift studies all agree that the 
M,/Mbuige relationship appears to be evolving with redshift. 
This evolution was seen to be independent of radio-loudness 



and quasar luminosity (see lDecarli et ai1l2010l . who studied 
both radio-loud and radio-quiet quasars, as well as inves- 
tigated possible biases). The results of these studies agree 
with one another and imply that the black holes at high 
redshifts are more massive for a given bulge mass than their 
local counterparts. 

Conversely, the op posite was found in 2 ~ 2 SMGs, 
d Alexander et al]|2008l ). where it was found that the SMGs 
lie below the local M,/Mbuige relationship. This illustrates 
that the selection of objects has significant effects as to 
where they will be placed with respect to the local relation- 
ship. SMGs are biased towards extreme star-formation rates, 
which indirectly translates to a bias towards relatively mas- 
sive galaxies. On the other hand, selection of quasars and ra- 
dio galaxies is biased towards massive black holes. It is there- 
fore perhaps unsurprising that these selected populations are 
biased towards different sides of the local M, /Mbuige rela- 
tionship. 

While our results have considerable uncertainty, they 
are consistent (wit hin the scatter) with what these other 
groups have found I Peng et a l. 2006b; McLure et al. I I2OO6I : 
iDecarh et alll2010l : Fargett. Dunlop fc McLu re 201^). For 
AMS12 to evolve to the local relationship, the bulge would 
have to grow ~ 3 times as much as the central black hole 
from 2 = 2.8 to 2 = 0. 

Given the amount of molecular gas implied by the CO 
observations, if AMS12 were to evolve secularly, it would, 
at most, only increase its bulge mass by ~ 15%. Mergers 
could add more gas for star formation while adding stel- 
lar and black hole mass. It is expected that massive galax- 
ies {M^, ^ 1 0^^ Mb), undergo > 1 mergers from 2 ^ 3 to 



les [Mi, fC l U Mb), undergo > 1 mergers trom 2 ^ 3 to 
present day jConsehce et al.ll2003l. [20071 : iBluck et al.ll2009l : 



S^^^^^^^^^^^^^^^^^^^^^^^^^, - - — . 

I2OIOI : iRobaina et aLlbOlOl ). However, it is not 
clear how mergers affect the black hole mass and whether 
it is possible to achieve the necessary growth of the bulge 
relative to the growth of the black hole. 

Alternatively, iDecarli et al.l (|2010l ) addressed the possi- 
bility that the remnants of the high redshift quasars they 
studied are high-mass outliers to the local relationship. 
These high redshift quasars are progenitors to present day 
massive ellipticals, and keeping their M,/Mbnige value to 
2 = 0, they become outliers, rather than evolve, to the local 
relationship. 

We must note that there was significant bias towards se- 
lecting a powerful quasar, giving a large M, , while demand- 
ing a faint 3.6 fj.m flux, limiting the host galaxy's luminos- 
ity (|Martfnez-Sansigre et "alll2005l . [2006^ 1. In addition, the 
search for CO in this object was initiated by selecting the 
brig htest MAMBO detection from the obscured quasar sam- 
ple jMartmez-Sansigre et~aLll2009l ). 

There is also a possibility that by biasing ourselves to- 
wards such a high Lboi, AMS12 is super-Eddington, in which 
case we would be overestimating the black hole mass. If 
AMS12 were accreting at super-Eddington rates, the black 
hole mass would be overestimated by the amount by which 
the bolometric luminosity exceeds the Eddington-limited lu- 
minosity. 
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6.7 Comparison to other galaxies 

AMS12 is the first unlensed, high redshift source detected in 
both [CI] lines. This eliminates any ambiguity on the effects 
of possible differential magnification from lensing. 

The characterisation of the gas and dust in AMS12 
is consistent with the observations of other high red- 
shift galaxies, including the value of log(LFiR/Lco), 
(kn own as the star format i on ef ficiency, i.e. Figure 8 
in ISolomon fc Vanden Bou^ l2005t ). where AMS12 lies 
within the scatter of the high redshift galaxies. The 
gas masses in these sources are similar to AMS12 (i.e. 
studie s of SMGs, ULIRGs and qu a sars at high red- 
shifts ISolomon fc Vanden Boud l2005l: ICreve et al.l 1 2005 : 



Riechers et al.l l2006l : ICoppin et al.l l2008l : lYan et al.l I2OIC : 
Lacv et al.ll201ll V 

A study of CO(l-O) in two z ~ 2.8 obscured quasar 
hosts found sim ilar gas and dust masses to AMS12 
l|Lacv et al.ll201ll ). The studies of obscured quasars and their 
hosts also indicate mature systems, with dust and gas masses 
low c ompared to the stellar mass estimates (e.g. iLacv et al.l 

However, comparing to other high redshift quasars 
which are strongly lensed, is difficult as is illustrated with 
F10214. The gas and dust properties of F10214 are very 
similar to AMS12, with approximately equivalent [CI] abun- 
dances, cooling rate, and line ratios in terms of both CO and 
[CI].^ 

lAo et all l|2008f) modelled the CO emission in F10214 
with LVG models, and found a similar dust temperature 
to AMS12 (80 K), and determined a range of gas kinetic 
temperatures from 45-80 K (note their Figure 7 has a similar 
shape to the high temperature region in Figure [5|, with Tex 
of [ CI] ^ 42 K. 

iRiechers et"ai1 l|2011ah detected CO(l-O) in F10214, and 
found there was no evidence for an extended, low-excitation 
gas component. For their analysis they assumed a con- 
stant magnification factor for the CO(l-O) (the magnifica- 
tion given by the higher- J CO lines), hence if there is dif- 
ferential magnification of the gas components, their results 
may be affected. 

Deane et al. (2012a, in prep.), have revised the lens 
model in this F10214 and have indeed found differential 
magnification on frequency and spatial scales. Deane et al. 
(2012b, in prep.), study resolved CO(l-O) in this object and 
find preferential magnification between individual channels 
and predict distortion of the CO SED. Thus, AMS12 which 
is unlensed, offers, so far, a unique opportunity to study the 
gas and dust in an obscured quasar host without the added 
complication of gravitational lensing. 



7 SUMMARY 

In this paper we have presented new observations of the 
obscured quasar AMS12 and, along with previous mm and 
submm observations, we have investigated the dust and gas 
properties of this object. 

The FIR dust observations are well fit by a single 
component graybody model with dust temperature Td = 
88 ± 8 K, emissivity index /3 = 0.6 ± 0.1, and Lfir = 
(3.2 ± 0.7) X 10" Lq, implying heating by the AGN. 



The CO SED was fit with LVG models, and we used 
the marginalised PDF of the dust temperature as the prior 
distribution of the gas kinetic temperature to constrain the 
parameters. This yielded the gas kinetic temperature of Tg 
— 90 K, and density n(Il2) = 10'^'^ cm~^, suggesting that 
SF is not the sole heating source. 

The atomic carbon fine structure lines [CI](^Pi — Po) 
and [CI](^P2 — ^ Pi), were observed and the [CI] excitation 
temperature was determined to be 43±10 K, which is signif- 
icantly lower than Tg, indicating either [CI] is not in LTE, 
or it is from a more extended, lower temperature gas com- 
ponent. 

The gas mass found from the CO(l-O) estimate, to be 
~ 4.9 X 10^" M0. The dynamical mass was calculated from 
the CO linewidth to be MdynSin^i = (2.7 ± 0.2) x 10^ Mq 
assuming r = 0.8 kpc, giving a limit to the host galaxy's 
inclination i 13°. 

The stellar mass in this object is estimated at M* = 
(3.2 ±0.3) X 10" Mq. It follows that; the gas and dust mass 
are only a fraction of the current stellar mass. The M./M* 
ratio is ^ 0.005, higher than in the local Universe. 

The system has already amassed the majority of its stel- 
lar mass and is host to a massive black hole, indicating a 
mature system. It is not clear how the system will evolve 
to the present-day Af,/Af* relation, or whether the extreme 
value is due to a selection bias. 
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